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ABSTRACT 

We present results of a 37 ks observation of CL 0016+ 16 with the XMM-Newton EPIC in- 
strument. Within 1.5 arcmin of the cluster centre we measure a gas temperature of kT = 
9.13^q 22 keV and an abundance of 0.22+QQg times the solar value (ler uncertainties). This 
significant improvement over previous measurements has allowed us to revise the estimate of 
the Hubble constant based on CL 0016+16 to 68 ± 8 kms -1 Mpc -1 (random error only), 
close to the value from the Hubble Space Telescope distance-scale project. The total gravi- 
tating mass within a radius of 248 kpc of the cluster centre is in good agreement with that 
found from gravitational lensing over the same region, supporting the assumption of isother- 
mal gas in hydrostatic equilibrium. The gas mass fraction of 0.13 ± 0.02 is in remarkable 
agreement with that given by cosmological parameters for the Universe as a whole, suggest- 
ing that CL 0016+16 is a fair sample of the matter content of the Universe. While there is no 
spectral or spatial evidence to suggest a cooling flow in CL 0016+ 16, we find an asymmet- 
rical central X-ray structure which may have a harder spectrum than the cluster as a whole, 
and be evidence for some merger activity, in addition to a previously reported asymmetry to 
the west of the cluster. The nearest companion cluster to CL 0016+16 is measured to have a 
gas temperature of kT = 3.8^0 3 keV and an abundance of 0.6^q'2 times the solar value (la 
uncertainties). We also present spectral data for the companion quasar whose spectrum was 
confused with CL 0016+16 in previous ASCA data. 
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1 INTRODUCTION 

The galaxy content of cluster CL 00 16+ 16 was described by 
Koo (1981) following its discovery by Kron (1980) in a photo- 
metric survey for faint galaxies. Its high redshift, z — 0.5455 
(Dressier & Gunn 1992), evident richness, high X-ray luminosity 
(L x , o.s-4.5 keV ~ 2 x 10 45 ergs s" 1 ; White, Silk & Henry 1981), 
and large Sunyaev-Zel'dovich effect (Birkinshaw, Gull & Moffet 
1981) led to extensive work on its properties as the prototypical 
example of a distant rich cluster. 

More recent X-ray and optical observations of the field of 
CL 0016+ 16 have found other nearby clusters of galaxies of lower 
mass but similar redshift (Hughes, Birkinshaw & Huchra 1995; 
Connolly et al. 1996; Hughes & Birkinshaw 1998a), suggesting that 
CL 0016+ 16 is the dominant member of a supercluster at z ~ 0.5. 
A quasar a few arcmin north of CL 0016+ 16 (Margon, Downes & 
Spinrad 1983) is also associated with this structure. 

The gas properties of CL 0016+ 16 have been the subject of 
particular study. Two reasons for this stand out: the intracluster gas 
can be used to study the cluster mass and baryonic mass fraction, 
and the combination of X-ray and Sunyaev-Zel'dovich data can be 



used to estimate the cluster's distance and hence the Hubble con- 
stant or other cosmological parameters. 

Grego et al. (2001) used X-ray data and the Sunyaev- 
Zel'dovich effect of CL 0016+ 16 to measure a gas mass fraction, 
/ g , that is close to the value seen in nearby clusters, implying that 
there is little evolution in f s between z ~ 0.5 and today. If clus- 
ter gas is the dominant baryonic component of CL 0016+ 16, and 
CL 0016+ 16 is a fair sample of the matter content of the Universe, 
/ g will be close to the baryonic mass fraction of the Universe, 
Sib / fi m • However, the temperature of the cluster gas is a necessary 
component of the calculation of both the gas and total masses, and 
improved measurements of the gas temperature would be expected 
to improve the reliability of cosmological conclusions drawn from 
studies of the mass components of the cluster. 

Good knowledge of the cluster temperature would also allow 
a comparison between the mass derived from the assumption that 
the cluster gas is in hydrostatic equilibrium and the mass measured 
by Smail et al. (1997) using a gravitational shear technique. This 
could set limits on the three-dimensional shape of CL 0016+ 16, 
which is a source of systematic error in the use of any cluster as a 
cosmological tracer. 
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Three independent Sunyaev-Zel'dovich effect measurements 
of CL 0016+ 16 (Hughes & Birkinshaw 1998b; Reese et al. 2000; 
Grainge et al. 2002) have been combined with X-ray data to mea- 
sure the distance of the cluster. All three lead to values for the Hub- 
ble constant of about 47 kms" 1 Mpc" 1 assuming qo = \ (this 
would be 57 kms" 1 Mpc" 1 for a cosmology with S7 m = 0.3, 
SI a = 0.7). The low value of the Hubble constant found by this 
technique might be attributed to projection or selection effects, cal- 
ibration uncertainties, unresolved substructure, or some other prob- 
lem, but it is clear that a substantial fraction of the random compo- 
nent of the error is due to uncertainty in the X-ray temperature. A 
deep X-ray observation with XMM-Newton, which can detect arriv- 
ing cluster photons of energy up to about 10 keV, and so can make 
a better measurement of the temperature of a hot cluster such as 
CL 0016+ 16 (which Hughes & Birkinshaw 1998b found to have 
kT ~ 7.55 keV in the rest frame), would be expected to resolve 
this problem. 

Thus we obtained a long XMM-Newton observation of 
CL 0016+ 16 with the aim of improving the measurement of the 
temperature and metal abundance of the cluster, and hence mak- 
ing it more useful as a cosmological tracer. This paper describes 
the data, their analysis, and the implications for the temperature 
and structure of CL 00 16+ 16 and one of the companion clusters 
in the supercluster. We also present the XMM-Newton spectrum 
of the companion quasar. We revise the estimate of the Hubble 
constant based on CL 0016+ 16 using our new temperature for 
the cluster, and compare an estimate of the total gravitating mass 
with that measured from gravitational lensing. We compare the gas 
mass fraction in CL 0016+ 16 with the value for the Universe as a 
whole given by S7b /S7 m . Finally, we point out structural features of 
CL 00 16+ 16 which require further study. 

We use a cosmology in which fi m = 0.3 and S7a =0.7, and 
we adopt H = 70 km s _1 Mpc" 1 outside the discussion of the 
use of CL 00 16+ 16 as a tracer of the Hubble flow. 



2 OBSERVATIONS 

XMM-Newton observed CL 00 16+ 16 on December 30th, 2000. In 
this paper we concentrate on data from the European Photon Imag- 
ing Camera (EPIC) which incorporates one pn (Sriider et al. 2001) 
and two MOS (MOS1 and MOS2; Turner et al. 2001) CCD-array 
cameras. The cameras were operated in Full Frame mode with the 
medium optical blocking filter. The data were provided to us in 
the form of two separate observations, of maximum livetime (for 
the MOS cameras) roughly 31.3 ks and 5.5 ks, separated by 110 
minutes. Analysis presented here uses the latest software available 
from the XMM-Newton project at http://xmm.vilspa.esa.es/ SAS 
V5.3.3. Calibrated event files were generated using the EMCHAIN 
and EPCHAIN scripts (described in detail on the project WWW site) 
and, after merging the data for each camera separately for the two 
observations, we have selected for analysis the good events with 
patterns to 12 from the MOS data and patterns to 4 from the pn 
data. 

The light curves extracted from the full-field, pattern=0, 10- 
12 keV events were examined for intervals of high particle back- 
ground, which we defined as those with > 0.1 ct s" 1 for MOS1 
and MOS2, and > 0.2 ct s" 1 for the pn. After removing the high- 
background intervals, the remaining exposure times were roughly 
35.9 ks, 36 ks, and 24.3 ks for MOS1, MOS2, and pn, respectively. 

Our analysis follows the prescription in Arnaud et al. (2002). 
This is superior to the simple use of local background in analysing 




Figure 1. Combined MOS1, MOS2 and pn data for the observation of 
CL 0016+16 in the 0.3-5 keV energy band, with a logarithmic intensity 
scale, a pixel size of 3 arcsec, and showing J2000 equatorial coordinates. 
Chip boundaries are visible. This image is vignetting corrected and with- 
out background subtraction. The bright source at a declination of roughly 
16°29'26", north of CL 0016+16, is the quasar E 0015+162 (Margon et 
al. 1983), and the poor companion cluster RX J0018. 3+1618 (Hughes et 
al. 1995) is visible to the SW. 



extended-source data because it treats separately the non-vignetted 
particle background and the vignetted cosmic X-ray background. 
First the EVIGWEIGHT task was run on both CL 0016+ 16 and 
blank-sky data so that after their subtraction the spectra and 
surface-brightness profiles (with local background) that are ex- 
tracted are corrected for vignetting. The vignetting-corrected 
MOS1, MOS2 and pn event files have been combined for display 
purposes. FigureQshows the 0.3-5 keV image of a wide field, cen- 
tred on CL 0016+ 16, from this merged event file. 

We use the data compiled by Lumb (2002), taken from deep- 
sky exposures cleaned of sources, to subtract from our data the par- 
ticle background, which is assumed to provide all the pattern=0 
counts between 10 and 12 keV. After applying our good-time- 
interval screening, the full-field, pattern=0, 10-12 keV count rates 
are the same in the CL 00 1 6+ 1 6 and blank-sky fields to within 1 %, 
7% and 16%, for the MOS 1, MOS2, and pn, respectively. The rele- 
vant scaling factor is applied to the vignetting-corrected blank-sky 
data before subtraction from the image data. 

The blank-sky data are not expected to be a good represen- 
tation of the background in the CL 0016+16 data at low energies. 
Here the cosmic X-ray background is the dominant contributor and, 
more importantly, at energies below about 1 keV the medium op- 
tical blocking filter used for our observations has a lower trans- 
parency than the thin filter used for the blank-sky data. This re- 
sults in an over- subtraction of background at the lowest energies. 
This is corrected for by using a blank-sky-subtracted source-free 
region of the CL 0016+ 16 image as a background for the blank- 
sky-subtracted on-source data. By selecting this background region 
to be an annulus of radii 4 and 11 arcmin, between position an- 
gles 90 and 120 degrees, and 220 and 310 degrees, we limited the 
number of source-contaminated exclusion regions to one circle of 
radius 30 arcsec. 

As data have been corrected for vignetting, spectral fitting uses 
on-axis response matrices and effective areas. We have used the 
version 6 response matrices made available by the project, and have 
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generated on-axis effective areas using the ARFGEN task. Bins be- 
tween energies of 0.4 and 10 keV are included in spectral fitting 
using the XSPEC software. 



3 EMISSION- WEIGHTED CLUSTER TEMPERATURE 

For good signal to noise, we have found the cluster spectral pa- 
rameters using the net counts (roughly 32,600 counts between 0.4 
and 10 keV) extracted from a circle of radius 1.5 arcmin. We fitted 
the data to a redshifted MEKAL model modified by local Galactic 
absorption. We let the relative normalization between the MOS1, 
MOS2 and pn data be free but all other parameters are in common. 
This results in MOS1 and MOS2 normalizations in good agree- 
ment, but the pn normalization is a few per cent lower, consistent 
with the evaluation of the EPIC flux calibration by Saxton (2002). 

We obtain a value of N H = (4.3 ± 0.6) x 10 20 cm -2 (95% 
confidence for one interesting parameter) for the Galactic hydrogen 
column density when it is a free parameter, in agreement with the 
value of 4.04 x 10 20 cm -2 found from 21 cm observations (Dickey 
& Lockman 1990). The value of the redshift fitted from the X-ray 
spectrum, z = 0.532±° °^, is slightly low as compared with pub- 
lished optical values of 0.5455 (Dressier & Gunn 1992) and 0.5481 
(Ellingson et al. 1998), suggesting that, although the EPIC energy 
scale is believed to be calibrated to better than 5 eV (Kirsch 2002), 
there may be a small additional systematic error in the CL 00 16+16 
observations. We have fixed the Galactic absorption to the 21 cm 
value but have allowed the cluster redshift to be free in finding the 
cluster parameters and uncertainties reported in this paper. 

The best fit gives kT = 9. 131^22 (±0.44) keV and an abi- 
dance of 0.22lQQ3(lQQg 5 ) times the solar value, where uncertain- 
ties are statistical and quoted as la (unbracketed), with 95% er- 
rors given in brackets, for one interesting parameter. The model is 
shown together with the data in Figure|2| and joint-confidence un- 
certainty contours are shown in Figure|3| With a reduced x 2 of 1 .09 
(X 2 = 859 for 786 de grees of freedom) the null hypothesis proba- 
bility is only 4%, but this is likely to result from the high statistical 
precision of the data coupled with small remaining systematic cali- 
bration uncertainties (Saxton 2002). An apparently random scatter 
in the residuals to the pn data is diminished if the pattern=0 events 
only are used for the pn, making the combined MOS and pn fit 
acceptable (reduced x 2 of 1.01: x 2 = 706 for 697 degrees of free- 
dom). However, as the best-fit parameter values are unchanged, we 
elect to include patterns 0-4 in the fits to improve statistical accu- 
racy. 

There is consistency between the parameter values fitted sep- 
arately for data from the MOS (kT — 9.3_tJ 7 , abundance = 
0.22±° ° l times the solar value, 2 = 0.534±g;^ ,'jVh = 4.4t?; u x 
10 20 cm" 2 ) andpn (kT = 8.6± ;^ , abundance = 0.22±g; O9 times 
the solar value, z = 0.525±^ 2 , = 4.3± u ;£ x 10 20 cm -2 ), all 
quoted with 95% uncertainties for one interesting parameter. 

The observed 2-10 keV flux within the circle of 1.5 ar- 
cmin radius is measured to be 1.43 x 10 -12 ergs cm" 2 s" 1 and 
1.39 x 10 -12 ergs cm" 2 s" 1 with the MOS and pn cameras, re- 
spectively. Between 0.4 and 10 keV the values are 2.24 x 10" 12 
ergs cm" 2 s" 1 and 2.18 x 10~ 12 ergs cm -2 s _1 , respectively. The 
average volume-weighted emission measures ( J n e n p dV / 'Ait D 2 ) 
are (1.22 ± 0.02) x 10 11 cm" 5 and (1.18 ± 0.02) x 10 11 cm" 5 , 
respectively (la uncertainties) 1 . 

1 In the XSPEC X-ray spectral-fitting software used here, and by most au- 
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Figure 2. XMM-Newton EPIC spectrum of CL0016+16 from counts 
within a circle of radius 1.5 arcmin. The upper spectrum is from the pn, 
and the lower spectrum is an overlay of MOS1 and MOS2 data points. The 
fit is to an isothermal model with kT = 9. 1 keV and an abundance of 0.22 
times the solar value. 
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Figure 3. CL 0016+16 spectral-parameter uncertainty contours using 
counts within a circle of radius 1.5 arcmin. Contours are lcr, 90% and 99% 
for two interesting parameters, with the normalizations and redshift allowed 
to vary. 

4 GAS DISTRIBUTION 

The intensity distribution of X-ray emission in the central part of 
CL 0016+16 (Fig.|4j shows an apparent filamentary structure in the 
shape of an inverted 'V. A true colour image, made using software 
from the Chandra project (the CIAO DMIMG2JPG task), with counts 
in 0.3-1 keV, 1-2 keV, and 2-8 keV corresponding to red, green 
and blue, respectively, suggests that the emission at the centre of 
the 'V may be hotter than in the surrounding part of the cluster. 
However the prong to the SE appears cooler than average. 

These spectral differences also appear when we use the 
IRAF/STSDAS ELLIPSE program to model the underlying large- 



thors, the normalization given by fitting thermal plasma models is 10 1 x 
(1 + z)' 2 times the volume-weighted emission measure. 
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Figure 4. Same as Figure 1, showing only the centre of CL 0016+16. Scale 
is linear with black corresponding to about 70 counts and the faintest pixels 
containing about 20 counts. An inverted V-shaped filamentary structure is 
apparent. This is also seen in the individual MOS and pn images, and is 
unrelated to chip orientation. 

scale structure of the cluster, and then subtract a scaled version of 
the model first from a 0.3-1.5 keV image and then a 2-8 keV im- 
age. The most prominent positive residuals in the soft image track 
the SE prong of the inverted ' V , whereas in the hard image they are 
at the cluster centre. The diffuse asymmetric emission to the west 
of the cluster, at about 1.5 arcmin from the centre (Fig. 1), appears 
somewhat softer than CL 00 16+ 16 as a whole. This western dif- 
fuse emission was previously discussed by Neumann & Bohringer 
(1997) using ROSAT data, and is interpreted as a merging subcom- 
ponent of the cluster that is associated with a peak in the weak- 
lensing mass map of Smail et al. (1995). 

We have fitted separately the counts within a source-centred 
circle of radius 0.6 arcmin and a source-centred annulus of radii 
0.6 and 1.5 arcmin using the method and model described in Sj3] 
The results (Fig. |5J indeed show that the inner region is hotter and 
of higher abundance than the outer one, at a combined confidence 
of 99% for two interesting parameters. The evidence for spatial and 
temperature structure in the inner regions may result from merger 
activity. A deeper observation would be required to probe this in 
detail. 

We have extracted the circularly-symmetric radial profile of 
background- subtracted 0.3-5 keV counts and fitted it to a spher- 
ical, isothermal, /3 model of surface brightness b x (0) oc (1 + 
[6»/6» c ] 2 )" 3,3+0 ' 5 , convolved with the point spread function (PSF). 
While the cluster is clearly ellipsoidal in appearance, and displays 
hotter regions within its core, the value of /3 and the average core 
radius of the model produced by fits of spherical functions are close 
to those obtained by making more complicated ellipsoidal fits (see 
Tables 1 and 2 in Hughes & Birkinshaw 1998b). If the fit results are 
to be used in the calculation of the distance of the cluster (as, for 
example, in S0, then values obtained from ellipsoidal fits could not 
be used directly to estimate Ho, since the projection is unknown for 
any individual cluster. However, Ho estimates obtained from spher- 
ical model fits for a population of clusters are almost unbiased, pro- 
vided that the clusters are selected in an orientation-unbiased fash- 
ion. Thus spherical modeling of clusters is sufficient for the purpose 
of building up a sample of distance measures from clusters. 

To construct the radial profile we have summed the counts 
from the three cameras. For the PSF we have used the on-axis pa- 
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Figure 5. Same as Figure 3, separating counts within a circle of radius 0.6 
arcmin (upper right contours) from those within an annulus of radii 0.6 and 
1.5 arcmin. 

rameter values for 1 keV from Ghizzardi's (2001, and private com- 
munication in 2002) fits of a /3 model to data from stellar sources. 
We have summed the MOS1, MOS2, and pn PSFs, weighted by 
each camera's contribution to the counts from CL 0016+ 16. The 
counts from a circle of radius 20 arcsec around the quasar 3.25 ar- 
cmin to the north of CL 0016+ 16 (Fig.[T} are excluded from the 
fits. 

The radial profile and best-fit model are shown in Figure |S| 
While the filamentary structure seen in Figure [4] leads to a poor 
fit within about 20 arcsec of the centre, this has negligible effect 
on the values found for j3 and 6 C , and the overall fit is accept- 
able (x 2 = 50.3 for 37 de grees of freedom). The values found 
are (3 — 0.697 ± 0.010, and 6 C = 36.6 ± 1.1 arcsec, where, since 
the parameters are strongly correlated, the errors quoted are one 
sigma for two interesting parameters. Nothing in the radial profile 
or spectral properties of CL 0016+ 16 is suggestive of a cooling 
flow. 



5 COMPANION CLUSTER, RX J0018.3+1618 

The companion cluster to the SW of CL 0016+ 16 (Fig.[TJ appears 
redder than CL 00 1 6+ 1 6 in a true colour image (made as in [J4j . We 
have found spectral-parameter values using on-source net counts 
from a circle of radius 1.5 arcmin (roughly 3,200 counts between 
0.4 and 10 keV) in a similar manner to the extraction and fitting 
for CL 00 16+ 16 itself. The fit to a MEKAL model is good, with 
X 2 = 109 for 101 degrees of freedom. The column density of 
A^h = 5.8l?.;3 x 10 20 cm" 2 (la errors) is consistent with the 
value from 21 cm measurements (4.09 x 10 20 cm -2 ; Dickey & 
Lockman 1990), and so we fix the value to the radio-derived mea- 
surement. FigureQshows the data and model whose best-fit param- 
eter values and uncertainties are as follows, where the latter are la 
(95%) confidence for one interesting parameter where all the other 
parameter values are allowed to vary: kT — 3.S^.Q'^(t.o'g), abun- 
dance = 0.Gt° i(t°oi) times solar, z = 0.472+°'°g (t°;°i>). The 
redshift is consistent (at 95% confidence, or \o if the systematic 
error of 0.0 1 6 suggested for CL 00 1 6+ 1 6 is applied) with the value 
of 0.5506 + 0.0012 found from optical measurements by Hughes et 
al. (1995) who were able only to place a lower limit on the cluster 
temperature (kT > 1.5 keV in the observed frame) using data from 
ROSAT. The precision with which the spectral parameters for the 
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Figure 6. 0.3-5 keV radial profile of CL0016+16 after subtraction of 
deep-sky and local backgrounds, as described in J2I The curve is the best 
fit circularly symmetric model convolved with the PSF: {3 = 0.697, 
9 C = 36.546 arcsec. The horizontal dotted line shows the contributions, 
taken into account in the fitting, of the model to the local background annu- 
lus. 
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Figure 7. XMM-Newton EPIC spectrum of the companion cluster from 
counts within a circle of radius 1 .5 arcmin. The upper spectrum is from the 
pn, and the lower spectrum is an overlay of MOS1 and MOS2 data points. 
The fit is to an isothermal model with kT = 3.8 keV and an abundance of 
0.6 times the solar value. 



companion cluster are measured is similar to that with which the 
parameters for CL 0016+ 16 were known prior to this work (Fu- 
ruzawa et al. 1998; Hughes & Birkinshaw 1998b). 

The net observed 2-10 keV and 0.4-10 keV fluxes are mea- 
sured to be 7.7 x 10" 14 ergs cm" 2 s" 1 and 1.6 x 10" 13 ergs 
cm -2 s _1 , respectively. The volume-weighted emission measure 
is (1.09 ±0.12) x 10 10 cm -5 . 

Because this cluster lies on a gap between CCDs in the pn 
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Figure 8. XMM-Newton EPIC spectrum of the companion quasar from 
counts within a circle of radius 40 arcsec. The upper spectrum is from the 
pn, and the lower spectrum is an overlay of MOS1 and MOS2 data points. 
The fit is to a power law of photon spectral index 2.38, with a narrow Fe 
fluorescence line of rest energy 6.4 keV. 



camera, as can be seen in Figure Q we have fitted the radial pro- 
file using only data from the MOS cameras. The parameter values 
found for a j3 model, convolved with the PSF corresponding to the 
off-axis position of the cluster, are of comparable value and pre- 
cision to the values found by Hughes et al. (1995) using ROSAT 
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6 COMPANION QUASAR, E 0015+162 

The companion radio-quiet quasar of redshift 0.554 (Margon et 
al. 1983) to the north of CL 0016+16 (Fig.[T) has a relatively steep 
spectrum, and we use the net counts between 0.3 and 10 keV within 
a source-centred circle of radius 40 arcsec for our fits. A fit to 
an absorbed power law is acceptable (\ 2 ~ H3 for 126 degrees 
of freedom), with the column density consistent with the Galactic 
value from 21 cm measurements. \ 2 ls reduced by 7.1 for one ex- 
tra parameter (significant at 99% confidence on an F test) when a 
narrow Fe fluorescence line, modeled as a Gaussian of fixed rest- 
frame energy 6.4 keV and a = 10 eV (small compared with the 
roughly 130 eV FWHM resolution of EPIC at the observed energy) 
is included in the model. If the line energy is allowed to be a free 
parameter, the best-fit value and la error is 6. lj^'f • Fixing the line 
energy at 6.4 keV, the following parameter values and uncertain- 
ties are found, where the latter are lcr (95%) confidence for one 
interesting parameter where all the other parameter values are al- 
lowed to vary: power-law photon index V = 2.38 ± 0.06(±0.13), 
iV H = 3.4 ± 1.0(±1.9) x 10 20 cm" 2 , equivalent width of the Fe 
line = 3081 2 23 (tin) eV. The spectrum is shown in Figure^ The 
observed 0.5-10 keV flux from the quasar is 1.61 ± 0.07 x 10" 13 
ergs cm" 2 s _1 (la error) and the luminosity is 2.15 x 10 44 ergs 

The AGN is roughly 8% the strength of CL 00 16+ 16 (in a 
cluster-centred circle of radius 4 arcmin) in terms of net counts 
between 0.4 and 10 keV (roughly 3,500 counts from the AGN). 
Reeves et al. (1997) fitted the ASCA spectrum of the AGN and 
found only an upper limit for narrow Fe fluorescence-line emission, 



© 0000 RAS, MNRAS 000, 000-000 



6 DM. Worrall & M. Birkinshaw 



but the ASCA spectrum would have been heavily contaminated with 
emission from CL 0016+ 16, and this would seem to account for 
the inconsistency between our results and the flat spectral index 
(1.99±0.09) and large flux (2.92 x 10" 12 ergs cm -2 s" 1 between 
0.5 and 10 keV) found by these authors. Our spectral index is con- 
sistent with the value of 2. 54 ±0.18 (lcr error) found for the quasar 
by Hughes & Birkinshaw (1998b) from a combined thermal and 
power-law model fitted jointly to ASCA and ROSAT PSPC data in 
order to characterize the emission both from CL 00 16+ 16 and the 
quasar. Their power-law normalization is also consistent with our 
better-determined value. XMM-Newton has improved the precision 
with which the spectral index is measured by a factor of three, and 
has provided the first evidence for an Fe fluorescence line in this 
radio-quiet quasar. 



al. 2002) of the effect. Using slightly different fi models for the 
cluster atmosphere, values A To = -1.20 ± 0.19, -1.24 ± 0.11, 
and —1.08 ± 0.11 mK were found in these investigations. We re- 
analysed the Hughes & Birkinshaw (1998b) data using the /3-model 
parameters based on our XMM-Newton data, and scaled the inter- 
ferometric data (which are relatively insensitive to the /9-model pa- 
rameters along degeneracy lines in the [/3, 9 C ] plane) to these same 
parameters. A weighted combination of all three independent val- 
ues leads to a central Sunyaev-Zel'dovich effect for CL 0016+16 
of AT = -1.26 ± 0.07 mK. Via equation Q, this implies a 
measurement of Daii c o- 

Combining our measurements of D\n^, and D\n e o, and us- 
ing the optical cluster redshift of 0.5481 (Ellingson et al. 1998), we 
find the angular diameter distance and central electron density of 
CL0016+16tobe 



7 HUBBLE CONSTANT AND GAS CONTENT OF 
CL 0016+16 

XMM-Newton has measured the temperature of CL0016+16's at- 
mosphere to ±2.5%. This is a considerable improvement on the 9% 
uncertainty obtained using ASCA and ROSAT (Hughes & Birkin- 
shaw 1998b), and so reduces one of the main sources of random 
error in the Hubble constant as deduced by comparing the X-ray 
and Sunyaev-Zel'dovich effect properties of the cluster. 

The XMM -Newton-derived emission measure within 1.5 ar- 
cmin of the centre, (1.20 ± 0.03) x 10 11 cm" 5 , from the weighted 
average of the MOS and pn values, implies an integrated cluster 
emission measure of 

J n e n p dV 

if the fraction of the emission contained within a 1.5-arcmin radius 
cylinder is extrapolated to the full cluster using the best-fit radial 
profile parameters. This emission measure is related to the angu- 
lar diameter distance of the cluster, Da, and the central electron 
density, n c o, by 



(1.75 ±0.05) x 10 11 cm" 5 



(1) 



EM = 



D A 9inl 0Fr(3/3-|) 



(2) 



J7(l±2) 4 4 T(3/3) 

if the XMM-Newton-deh\ed [3 model is used to describe the overall 
structure of the cluster, and so the X-ray data provide a measure- 
ment of Da^-oo- The symbol 77 in equation (2) is the electron/proton 
ratio, ~ 1.17. The emission measure Q is less than half the value 
given by Hughes & Birkinshaw (1998b), because of an error in the 
way that plasma thermal emission normalizations were reported by 
XSPEC. Since the emission measure was not used in the calcu- 
lation of the cluster distance in the Hughes & Birkinshaw paper, 
however, the conclusions about the cluster distance there are un- 
changed by this error. 

The (3 model can also be used to calculate the central Sunyaev- 
Zel'dovich effect for the cluster, 



ATo 



-2 T rad Or 



(^) 



n c o Da dcy/n 



r(|/3- 



(3) 



where T ra d is the thermodynamic temperature of the microwave 
background radiation and or is the Thompson scattering cross- 
section. 

The central Sunyaev-Zel'dovich effect of CL 00 16+ 16 is 
found by fitting a one-dimensional profile (Hughes & Birkinshaw 
1998b) or a two-dimensional map (Reese et al. 2000; Grainge et 



(4) 
(5) 

= 0.3, Qa = 0.7, this corre- 
± 8 kms" 1 Mpc" 1 (random 



D A = 1.36 ±0.15 Gpc, and 
n c0 = (8.8 ±0.5) x 10~ 3 cm" 

In a ACDM cosmology with fl n 
sponds to a Hubble constant of 1 
error only), which is about 45% larger than the values estimated by 
Hughes & Birkinshaw (1998b), Reese et al. (2000), and Grainge 
et al. (2002). The change in Hq arises about equally from the in- 
creased best-fit gas temperature found in the present XMM-Newton 
study and the change from a flat matter-dominated cosmology to a 
flat ACDM cosmology. 

In addition to the statistical error in equation 0, the distance 
is subject to systematic errors, the largest of which are 5% cali- 
bration uncertainties for the Sunyaev-Zel'dovich effect and X-ray 
normalization, a 20% uncertainty from possible projection effects, 
and a 15% uncertainty from substructure within the cluster (Hughes 
& Birkinshaw 1998b). Overall, therefore, we deduce an angular di- 
ameter distance for CL0016+16 of 1.36 ± 0.15(statistical) ± 
0.34 (systematic) Gpc, which implies a Hubble constant of 68 ± 
8(statistical) ± 18(systematic) kms" 1 Mpc" 1 . This is consis- 
tent with the value of the Hubble constant found in the HST 
distance-scale project (Mould et al. 2000), but the systematic er- 
rors are too large for a useful independent measurement of Q A to 
result from our distance for CL 0016+ 16. 

It is interesting to use the X-ray data for CL 0016+ 16 to es- 
timate the mass of the cluster within the region well-probed by 
XMM-Newton. For an isothermal gas in hydrostatic equilibrium in 
a spherical cluster, the total gravitating mass within a cylinder of 
radius r about the cluster centre is 



M t , 



3ivf3 
t{r) = -r- 



kT 



2 (i m p G r c 



2 / 2 
(l + ^ 



(6) 



For r — 248 kpc, this corresponds to a mass of (2.0 ±0.1) x 
10 14 M Q . This value is consistent with the mass of (2.7 ± 0.9) x 
10 14 Mq measured from the gravitational lensing effect of the clus- 
ter within the same cylindrical radius (Smail et al. 1997; value con- 
verted to our adopted cosmology). 

We can compare this total gravitating mass with the gas mass 
within the same cylindrical radius, based on the central electron 
density that we measure. The result, Af gas (248 kpc) = (2.6 ± 
0.2) x 10 13 M© indicates that X-ray emitting gas contributes a 
fraction (0.13 ± 0.02) of the total mass of the cluster in this 248- 
kpc radius cylinder (for Hq = 70 kms" 1 Mpc" 1 ), or (0.13 ± 
0.04) independent of cosmology if we use the cluster distance in 
equation and the electron density in equation 

The measured gas mass fraction is consistent with the value 
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0-14_o'q4 within a radius of 65 arcsec (415 kpc) which is obtained 
by correcting the result of Grego et al. (2001) to the cosmology 
that we use and cluster temperature that we measure. Our tem- 
perature makes CL0016+16's gas mass fraction more consistent 
with the average in the Grego et al. sample, and with the aver- 
ages for the Myers et al. (1997) and Mohr et al. (1999) samples. 
Since these samples derived the gas mass fraction using different 
analyses of the available X-ray and/or Sunyaev-Zel'dovich effect 
data, their agreement about the gas mass fraction and its lack of 
change with redshift suggests that the assumptions used, notably 
the lack of small-scale density or thermal structure in the atmo- 
spheres, are not seriously in error. Since the baryonic matter con- 
tent of the X-ray gas dominates that of galaxies in clusters, and 
fib /fi m = 0.12 ± 0.02 (Turner 2002) for the Universe as a whole, 
it appears that our 248-kpc radius cylinder through CL 0016+ 16, 
and clusters of galaxies in general, are fair samples of the matter 
content of the Universe. Studies of how the gas mass fraction in 
clusters varies with redshift are therefore unlikely to provide use- 
ful information on cluster evolution, at least at the current level of 
accuracy of the measurements. 



8 SUMMARY 

With the EPIC camera of XMM-Newton we have measured 
the temperature and distribution of the intra-cluster medium of 
CL 0016+ 16 to unprecedented precision. The emission-weighted 
temperature within 1.5 arcmin of the cluster centre is kT — 
9.13±J^ J| keV and the abundance is 0.22^°;°* (la uncertainties). 
The higher best-fit temperature than found in previous work us- 
ing ASCA (Furuzawa et al. 1998; Hughes & Birkinshaw 1998b) 
affects the cluster's position on the temperature-luminosity plane, 
now moving it into good agreement with the L x — T relation for 
nearby clusters (see figure 7 of Schindler 1999). 

Although there is some ellipticity in the plane of the sky, as 
previously reported based on ROSAT data (Neumann & Bohringer 
1997; Hughes & Birkinshaw 1998), we find an acceptable fit of 
the overall envelope out to 5 arcmin (1.9 Mpc) to a circularly- 
symmetric (3 model, representative of gas in hydrostatic equilib- 
rium, with (3 = 0.697 ± 0.010, and core radius 6 C = 36.6 ± 1.1 
arcsec (where errors are la for two interesting parameters since /3 
and 9 C are strongly correlated). 

While there is no spectral or spatial evidence to suggest a 
cooling flow in CL 0016+ 16, an asymmetric enhancement to the 
west of the cluster is apparent in the XMM-Newton image. Such 
an asymmetry was previously reported by Neumann & Bohringer 
(1997) and interpreted as a merging subcomponent of the cluster. 
The XMM-Newton data also show a more central asymmetric X- 
ray structure which may have a harder spectrum than the cluster as 
a whole, and be evidence of merger activity. A deeper X-ray expo- 
sure would be required to probe this in more detail. 

The improved spectral and spatial data provided by XMM- 
Newton have allowed us to reduce the random error on the angu- 
lar diameter distance of CL 00 16+ 16 to 11%. The resulting es- 
timate of the value of the Hubble constant, 68 ± 8 (random) ± 
18(systematic) kms" 1 Mpc -1 , is in good agreement with the 
value 71 + 6 kms -1 Mpc -1 from the HST distance-scale project 
(Mould et al. 2001). The largest remaining systematic component 
of the error in our estimate of the Hubble constant arises from pro- 
jection and substructure effects, and could be reduced by a fac- 
tor ~ 3 by applying the technique to a sample of 10 — 20 clus- 



ters with equally-good X-ray and Sunyaev-Zel'dovich effect data 
as CL 0016+16. 

The total gravitating mass within a radius of 248 kpc of the 
cluster centre is in good agreement with that found from gravita- 
tional lensing over the same region (Smail et al. 1997), which sup- 
ports the assumption that the gas in CL 00 16+ 16 can be treated as 
isothermal and in hydrostatic equilibrium. The measured gas mass 
fraction of 0.13 ± 0.02 is in remarkable agreement with that given 
by the cosmological parameters for the Universe as a whole, sug- 
gesting that CL 00 16+ 16 is a fair sample of the matter content of 
the Universe. 

The cluster CL 0016+ 16 is believed to be the dominant mem- 
ber of a supercluster at z ~ 0.5. The nearest companion clus- 
ter and relatively bright companion quasar E 0015+162 (which 
is of an X-ray luminosity more typical of a local Seyfert galaxy) 
are both well measured in the XMM-Newton EPIC data. In pre- 
vious ASCA data the emission from the quasar is confused with 
that from CL 0016+ 16 and so values reported for the quasar based 
on that observation are unreliable. We find that the spectrum of 
E 0015+162 fits well a power law absorbed only by the column 
density of gas through our Galaxy, with a relatively steep power- 
law photon index of T — 2.38 ± 0.06. There is evidence for the 
detection of an Fe fluorescence line. 

XMM-Newton finds that the spectrum of the nearest compan- 
ion cluster, RX J0018.3+1618, fits a thermal model with kT = 
3.8lg'3 keV and abundance = 0.61q 2 times solar, These uncer- 
tainties for RX J0018. 3+1618 are comparable to those of earlier 
measurements of CL 00 16+ 16. Deeper and broader X-ray obser- 
vations of the supercluster may provide important information on a 
massive structure at an early phase of its evolution. 
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